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The effect of Rashba spin-orbital interaction upon supercurrent in S-2DEG-S proximity junctions 
is investigated in the clean limit. Generalization of Beenakker's formula for Andreev levels to the 
case of spin-orbital scattering is presented. Spin-orbit - induced splitting of Andreev bound-states 
is predicted for a junction of infinite width, and with non- vanishing normal backscattering at S/N 
interfaces. Semiclassical average of the Josephson current is however insensitive to the Rashba 
coupling as long as electron-electron interaction in 2DEG is neglected. 
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I. INTRODUCTION 

Josephson junctions of two superconductors via two- 
dimensional electron gas, (most usually implemented 
in the Nb/InAs/Nb structures) were actively studied 
both experimentally and theoretically, cf. e.giSiMi^. 
Generic feature of all these devices is strong reduction of 
the experimentally measured product IqRn with respect 
to the theoretical predictions. In particular, this dis- 
crepancy is known for short junctions with high-quality 
S/N interfaces, demonstrated by measurement of a non- 
sinusoidal current-phase relation^. At the temperatures 
much below 7^, the parameter IqRn ~ 0.22mV was mea- 
sured in Ref. |a, to be compared with the Niobium super- 
conductive gap A « 1.5meV. Thus it seems natural to 
look for some effects which were not taken into account 
in the existing theory, cf. e.g.^'^, but could be responsible 
for such a drastic suppression of the critical current. 

An obvious candidate to be explored is the Rashba 
spin-orbital interaction^ Hji = Q;[f7 x p] • n known to 
exist in the 2DEG structures due to up-down asymmetry 
of the quantum well (here n is the unit vector normal 
to the plane of the 2DEG) . In the InAs heterostructures 
this term is especially large (cf. Ref . . leading to the 
band splitting An = 2apF ~ hmeV , i.e. considerably 
larger than the Niobium superconductive gap. Therefore 
it seems natural that the account of the Rashba term 
might be important for analysis of the Josephson current 
in these devices. We also note in this relation the paper—, 
where it was shown that persistent currents in mesoscopic 
metal rings should be modified strongly by SO coupling 
- which seems to point out the existence of similar effect 
upon the Josephson current. 

However, it is frequently assumed that the spin-orbital 
interaction cannot influence proximity effect in supercon- 
ductive structures, since it respects time- reversal invari- 
ance. This argument is not valid, however, when the crit- 
ical Josephson current is considered, since the presence of 
a current already breaks down the time-reversal symme- 
try. More detailed arguments seem to come from recent 
papersiSii^, where the influence of both Rashba coupling 
and Zeeman magnetic field upon the critical current of 
S-N-S junctions was considered. In both these papers it 
was found that in the absence of the Zeeman term the 



Rashba interaction (if treated within the simplest model 
of equal Fermi velocities on both chiral branches) totally 
cancels out from the equations for Andreev levels. We 
will show, however, that this cancellation is not generic; 
rather it is due to different simplifications used in the pa- 
pers mentioned: a model of completely transparent S/N 
interfaces was employed in Ref. |l2L and a purely one- 
dimensional model used in Ref. Il^ 

It will be shown below that in the general case, when 
some normal backscattering at an arbitrary angle of in- 
cidence occurs at the S-N boundaries, the spin-orbital 
coupling does affect the energies of the Andreev levels 
and the supercurrent they carry on. We will show that 
the effect of SO interaction can be understood as being 
due to modification of the transmission channels defined 
by the scattering matrix 5, which describes the junction 
properties in the normal state. For a model junction 
with infinite length (or periodic boundary conditions) in 
the direction transverse to the supercurrent, spin-orbit 
splitting of transmission eigenvalues is found, which re- 
sults in splitting of each Andreev level into pair of spin- 
polarized levels, with phase-dependent energy difference 
5E{x)- Note, that 5E{Q) = 0, in agreement with the 
time-reversal invariance restored in the absence of the 
phase bias. The idea that the Andreev levels can be spin- 
splitted due to the SO coupling was proposed in Ref. ^3 
for a narrow (few-channel) junction. The SO effect we 
discuss in the present paper is different from the one pre- 
sented in Ref. [3 

In this paper we consider the simplest two-dimensional 
model of baUistic S-2DEG-S junction (cf. e.g. Ref. 8) of 
infinite width in the lateral direction transverse to the 
current flow, see. Fig. 1. We neglect possible poten- 
tial barriers at the S/N interfaces, assuming that the 
normal backscattering is due to Fermi-velocity mismatch 
only, and consider ballistic electron propagation along 
2D structure between superconductive terminals. In the 
Sec. II of the paper we show that in the short junction 
limit (junction length i <C = Hvp/ A, where vp is the 
Fermi-velocity of 2DEG) the positions of the Andreev 
levels can be expressed via the transmission eigenvalues 
T of the full scattering matrix S in precisely the same 
way as was found by Beenakker for junctions with 
spin-independent scattering. Then in Sec. Ill we present 



2 



-L/2 



Rashba 
2D 



L/2 



FIG. 1: Two-dimensional model of a superconductor/Rashba 
2DEG/superconductor Josephson junction infinite in direc- 
tion perpendicular to the current (along y axis). The Rashba 
2DEG region has thickness L\ m/mn is the efi'ective mass and 
Va/vn is the Fermi velocity in the S/2DEG; ip is the angle be- 
tween the velocity direction of a quasiparticle and the x axis 
in the 2DEG region; n is a unit vector normal to the plane of 
the 2DEG. 



calculations of the scattering matrix S for the simplest 
two-dimensional model of ballistic S-2DEG-S junction 
(cf. e. g. Ref. 1^ of infinite width in the lateral di- 
rection transverse to the current flow, see. Fig. 1. We 
demonstrate explicitly spin-splitting of the transmission 
probabilities T±(py) for the transmission channels char- 
acterized by the momentum component py. We show 
then, that the distribution function for the transmission 
probabilities 7^(T) coincides with the one discussed by 
Melsen and Beenakker^* in the absence of spin-orbital 
coupling. In the Sec. IV we derive expression for the 
Josephson current of a short junction and demonstrate 
that average current is insensitive to the Rashba cou- 
pling. In Sec.V we go beyond the short-junction limit: 
we derive an equation for spin-splitted Andreev levels for 
the junction with arbitrary L/^q ratio and demonstrate 
that their contribution to the average (semiclassical) su- 
percurrent is insensitive to the SO coupling. Section VI is 
devoted to the discussion of applicability of our results to 
S-2DEG-S junctions of finite width and possible ways to 
detect spin-splitted Andreev levels. Finally, in Sec. VII 
we present our conclusions and discuss open problems; 
in particular, it is proposed that the account of electron- 
electron interaction together with Rashba coupling might 
be able to explain Josephson current. 



II. SPECTRUM OF ANDREEV LEVELS 

The excitation spectrum consists of the positive eigen- 
values of the Bogolyubov-de Gennes (BdG) equation: 
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where (C/*)"^ — g''°'[U{v, ^))* g^p, a, (3 are spinor in- 
dices, g = i&y is a metric tensor in spin space; ^ = 

2 

^ — Ep is the kinetic energy of a quasiparticle (energies 
are measured relative to the Fermi energy); 
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In our model in the normal region (mapped on Fig. 1 as 
"Rashba 2DEG") the operator U = a[a x p] ■ n is the 
spin-orbit interaction, which preserves the time-reversal 
invariance. In the superconductors the Rashba term is 
absent, U — 0. The superconducting gap A is assumed 
to be a step-like function: equal to zero in the normal 
region, and its' modulus |A| constant and equal in both 
superconductors. 

Equation which relates the excitation spectrum of the 
Josephson junction to the scattering matrix in the normal 
state S, was derived in Ref. 
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where 



rhe = irA, reh = "fr^, rA 
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rhe is the Andreev reflection matrix for e — s- /i scatter- 
ing in the space of channels incident (reflected) on the 
left and right NS boundary, ±x/2 are the phases of the 
left (right) superconductor, S^^h) is the electron (hole) 
scattering matrix of the normal state. 

When no spin-dependent scattering is present, normal 
scattering matrix Se is trivial in spin space, i. e. pro- 
portional to the unit matrix (Tq. Further on, within the 
short junction limit L <^ ^q, scattering matrices Se,h do 
not depend on energy, moreover Sh — S*. Therefore 
Eq.(|3Jl can be transformed to an explicit solution^^ for 

spin-degenerate Andreev levels, ej — ±^J I — Ij sin^ ^, 
where Tj is the j-th eigenvalue of the transmission prob- 
ability matrix T^T (eigenvectors of this matrix define 
scattering channels). Below we show that solution of the 
same kind can be obtained when the spin-orbital scatter- 
ing is present. 

In the presence of spin-orbital interaction the scatter- 
ing matrix in ^ becomes spin-dependent but still obeys 
time-reversal invariance. It makes possible to general- 
ize Beenakker's derivation for Andreev levels in short 
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using the following set of relations for the S- 



Sh{e,Py) = f Sl{~€,~py)g, 



(5) 



where the superscript ^ means full matrix transposi- 
tion. The first relation in Q is just a unitarity condi- 
tion, the second one follows from the time-reversal in- 
variance (we used here a transformation of the wave- 
functions upon time-reversal, ip*~^(py) = gip*{—py)). Fi- 
nally, the third relation in (jSJ is due to a special symme- 
try of the Bogolyubov-de Gennes equations: iphi^,Py) — 
g'^ipl{—€,~py). It is important to note the sign change 
of the parameter py in the second and the third relations 
above: when all scattering states are characterized by a 
conserved momentum (pj,), the time-conjugation opera- 
tion contains complex conjugation and py —py inver- 
sion, since time reversal of the scattering matrix should 
change the sign of px only, while keeping py intact. In 
other terms, the additionally —^ —py operation is needed 
due to the use of scattering channels characterized by 
complex eigenfunctions cx e*^"^. Usually in calculations 
of this kind a real basis of transmission channels is used, 
in which case such an additional operation would be ab- 
sent. 

Making use of relations one can transform Eq. Q 
into the form: 
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For a short contact L <C foi we neglect the energy- 
dependence of the scattering matrix in Eq. and ob- 
tain a second order equation with respect to e^, which 
results in the following solution: 



<^s,ri{Py) = 77AWI - Ts(p; 
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where ry = ± and 71 (py) are transmission probabilities 
- eigenvalues of the matrix T^^T, depending on the spin 
index s = ± and the conserved momentum py. In gen- 
eral, T+(py) 7^ T^{py), thus four non-degenerate Andreev 
levels correspond to each py value, as shown in Fig. 3 be- 
low. Note however, that full family of Andreev levels still 
contains pair-wise degeneracy within our model. Namely, 
degeneracy exists between states with py = ±|py|. Be- 
low we consider a specific example of scattering problem 
relevant to S-2DEG-S structures, and calculate Ts{py) 
eigenvalues. 



due to Fermi velocity mismatch only, Vs ^ Vn (here 
Vs and Vn are the Fermi velocities in the superconduc- 
tive metal and in the 2DEG correspondingly). Addi- 
tional source of reflection due to an effective potential 
barrier at the interface (cf. e. g. Ref. @) can be 
present, but docs not affect our results qualitatively. 
Since the effective mass rn„ in the 2DEG differs strongly 
from the effective mass m in the metallic superconduc- 
tor (typically, m„/m w 0.03 for 2D structures with 
InAs), the difference of these masses should be taken 
into account explicitly. Our first goal now is to find 
refiection/transmission amplitudes on single S/N inter- 
faces (for the normal state of the superconductive metal 
S). We will follow Ref. 0, and use continuity equations 
which follow from the Schrodinger equation with space- 
dependent mass m(x) and spin-orbital parameter a{x): 
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where F denote F{x = 0) - F{x ^^^-0) for 

the left interface (cf. Fig. 1) and similarly for the right 
interface located at x = L/2. Further on, Pp = mvs and 
Pf = runVn are the Fermi-momenta in the S metal and 
in the 2DEG correspondingly; usually ^ Pp, whereas 
Vs and Vn are of the same order of magnitude. Below we 
will assume that the parameter a/w„ <^ 1 measuring the 
relative strength of the Rashba interaction is small in 
comparison with the Fermi-velocity mismatch, i.e. a <C 
\vs — Vn\- Under this condition, reflection amplitudes 
at each of the S/N boundaries are determined by the 
ratio Vn/vs only. Then the amplitudes of reflection and 
transmission are trivial in the spin space, e.g. 'r'"^ — 
S"^T*i, and so on. For an incident wave incoming from 
a; = — 00, the amplitudes of reflection and transmission 
on the left (1) interface, are: 
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where w — Vnxj'Vsx is the ratio of the x-components of 

the electron velocities, with Vnx — "Vn cos and Vsx — 
1/2 

\p1 ~ (i^)^^n sin^ (/?] « Vs. Here ip is the angle be- 
tween the velocity direction and the x axis in the 2DEG; 
note that Vgx is very close to Vg for any angle ip, since 
(mn/rn)^ <C 1. Other reflection/transmission amplitudes 
are determined as follows: 



r 1= r 2 



r 2 = r I, t 2 = t I, 

<_ 2u; 
T^i, ii = t2^-—- 10 

l + w 



III. S-MATRIX AND TRANSMISSION 
EIGENVALUES 

We are interested in the study of speciflc spin-orbital 
effects and thus will consider the simplest model of S /N 
boundaries, assuming that normal electron reflection is 



The total scattering matrix S of the S/Rashba 
2DEG/S junction in the normal state, formed out of 
the "single boundary" amplitudes, Eq. reads (sim- 
ilar equations can be written for T2 and i?2): 

fi = Ta^'- [1 - Yi{S')-^T^2S'-] Ti, 
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(11) 



where R and T are the reflection and the transmission 
blocks of the scattering matrix: 



Rl T2 



(12) 



the index 1 in the amplitudes R, T in Eq. Hll() means 
that the equations are written for the case of an electron 
propagating from left to right. Matrices in the spin space 
describe spin rotation during the electron propaga- 
tion across 2DEG region with the Rashba coupling be- 
tween the two S/N boundaries. Explicit form of these 
matrices can be obtained by transformation of the plane- 
wave eigenmodes with definite chirality to the spin basis 
with definite Sy projection: 

_ j(.Qgy^ _ i sin A sin tpcTx + i sin A cos ipa z] (13) 
(o )^^ — e*^ [cos A — i sin A sin ifi&x — i sin A cos tpaz] ■ 

Here f(e) — k{e)L, with fc(e) = k + me/k and k — 
p F cos Lp, is the main semiclassical phase and A — 
m„aL/ cos (/? is the additional phase due to spin rota- 
tion by Rashba coupling. Within our approximation 
ct/vn ^ 1, the whole effect of the Rashba coupling is 
contained in the phase A which is not small if length L 
of the junction is comparable or larger than the spin- 
rotation length Lo = fi/mnCt. 

We define, for further convenience, a new parameter 
X = \ogj^, where w is defined after Eq. @. Then, 
combining Eqs. (|ll|l and (|13|l . we obtain a transmission 
matrix in the form: 



Ti = To + TifTj; -I- TsCTz, 
T2 = To + TiCTj. - T-i&z, 



(14) 



with 



Tq = t sinh(a; — i^) cosA, 

Ti = —it cosh(x — i^) sinA siniy9, 

T3 — it sinh(a; — i^) sinA cos ip, 

where we denoted 

sinha; 

t 



sinh (x — i^) + sin A sin ip 



(15) 
(16) 



The reflection matrix R has the form: 

i?l = Rq + Rl(Jx + R2(^yi 
R2 = Ro + RlO'x — R2'^y, 



(17) 



with 

Ro — t 



cotli X sin A sin ip 



. sin^ 
sinha; 



sinh(x — i^) 




FIG. 2: The spin-splitted transmission eigenvalues Ts, s = 
±1, as functions of the angle of propagation ip of the quasi- 
particles inside the 2DEG, plotted for a realistic S-2DEG-S 
junction with parameters Vs = 7 ■ lO^cm/s, v„ — 5 ■ lO^cm/s, 
m = nie, ra„ = 0.035me, mna/h = 5 • lO^'cm"^, L = 190nm. 
For these parameters and for value of the superconducting 
gap A = 1.4meV: (1) the length of the contact L is shorter 
than the coherence length, ^0 = hvs/A — 330nm; (2) the 
Rashba velocity is much smaller than the Fermi velocity in the 
2DEG, a/vn « 0.03; (3) the system is within the semiclassical 
limit, ppL/h = m„VnL/h ~ 30; (4) the spin-orbital splitting 
2ctpF ~ 3.3meV is larger than the superconducting gap A; 
(5) the S/N interfaces are almost transparent {vs/vn ~ 1-4), 
which allows a large experimental value of the critical current. 



Rl — —t sin2A simp, 

i , . 
i?2 = - i sin A sin2v3. 



(18) 



Now we use Eqs. H14I I15|l to obtain the transmission 
probabilities, as the eigenvalues of the matrix T = T'^T: 



sinh X 



sinh^a; + sin^(C± §)■ 



where the phase (3 defined via 



cos (3—1 — 2 sin ip> sin A 



(19) 



(20) 



is due to the Rashba interaction. Eq. H19|l demonstrates 
explicitly spin-splitting of the transmission eigenvalues 
T±. Note that (3 = and the splitting is absent for 
trajectories with tp = 0, which is the case of a purely 
1-dimensional (single-channel) contacti^. In the absence 
of normal reflection, i.e. at a; — > 00, all transmission 
eigenvalues are equal to unity and the spin-orbital effects 
disappear as well^^. 

Spin-orbital effect upon 7± reduces, according to 
Eq. H19|l . to the shift of the main semiclassical phase, 
^ — > ^ ± P/2, i n ag reement with the result presented in 
Eq. (1) of Ref. [ij. An example of the dependence of 
transmission eigenvalues T± as functions of the incidence 
angle is shown in Fig. 2. An important point to mention 
is that this dependence is even with respect to — > — 
reflection, cf. Ea. (|20|l . This symmetry reflects a "trace" 
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of the Kramers degeneracy which is known to exist for 
the transmission eigenvalues defined within the real basis 
of scattering states (note that the proof of Kramers de- 
generacy of transmission eigenvalues is by far more com- 
plicated!'' than original Kramers theorem for the degen- 
eracy of energy levels). We characterize scattering states 
by complex travelling waves e^^"^, which leads to break- 
down of time invariance. Therefore Kramers theorem is 
not applicable to our model and spin-splitting demon- 
strated in Ea. l(T^ may occur—. 

It follows from Eq. (|19|l that a semiclassical average of 
any physical quantity which can be expressed as a sum 
of terms containing individual T variables (i.e. does not 
contain cross-terms like 71)^71 ) , does not depend upon the 
SO coupling. Indeed, calculation of any average quantity 
in our model involves integration over momentum com- 
ponent py parallel to the interfaces (or on propagation 
angle Lp^ defined as Py = ppsmip). The integrand, as 
function of ip, contains fast oscillations with characteris- 
tic scale 1/pfL, and relatively slow dependence on cos (p. 
It is convenient first to average over fast oscillations by 
going to probability distribution of transmission eigen- 
values defined as 



(21) 



Clearly, the presence of ±/3 phase-shift does not alter the 
form of the probability distribution, which is of the same 
form as considered, e.g. in Rcf. 18 and independent on 
the spin-orbit coupling: 



tanh : 



2TV1~T - tanh^ 



(22) 



Now we consider, as the simplest example, calculation 
of the average conductivity of a junction in the normal 
state. It can be written as 



G G: 



cosip^ f TV^{T)dT. 

J-7T/2 J 



(23) 



Universality of the distribution fmiction Vip{T) leads to 
independence of the average conductance G, as well as of 
other quantities which can be expressed via this distribu- 
tion functions, upon the spin-orbit phase A (remember 
that we neglected weak effects of the order oi a/vn ^ 1 
). Note once again, that the above simple considera- 
tions could not be applied to calculation of any quantity 
which is not additive as a function of different transmis- 
sion channels, i.e. which contains products of different 
transmission eigenvalues. 



IV. JOSEPHSON CURRENT 

Equation for the Andreev levels (O together with 
Eq. H19|) for transmission eigenvalues constitute the cen- 
tral result of the paper. Now one can calculate Josephson 
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FIG. 3: The four spin-splitted Andreev levels its, s — ±1, as 
a function of the superconducting phase difference x, plotted 
for a value of the angle of propagation = tt/5, and for a re- 
alistic S-2DEG-S junctions with parameters Vs = 7 ■ lO^cm/s, 
Vn — 5 ■ lO^cm/s, a ~ 0.2 ■ lO^cm/s, m — nie, rUn — 0.035me, 
L = 190nm. 
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Ts{Py) 

(x) 



tanh 



(X) 



2T 



(24) 

Eq. H24|l is applicable for the temperature-dependent 
Josepshson current in the short-junction limit. In the 
semiclassical limit ppL oo the average Josephson cur- 
rent can be calculated with the use of the distribution 
function V(p{T) given by Eq. (|22|l as follows: 



iix) = 



eA 

2h 



If COS Lp 



7r/2 



X (25) 



V^{T)dT 



Tsinx 



l-Tsin^l 



: tanh 



Ajl -Tsin^f 



2T 



Eq. (|25|l demonstrates independence of the average 
Josephson current on the spin-orbital coupling. Such an 
average current is a meaningful characteristic of a junc- 
tion with both lateral sizes much longer than the Fermi 
wavelength, L,Ly ^ h/pF- 

Oscillations of Ic as a function of the electron density, 
was discussed theoretically in Ref. within the model 
very similar to the present one (but without Rashba cou- 
pling) . It was argued that oscillations should appear due 
to the presence of normal resonances in a double-barrier 
structure, like the ones described in our Eq. (|19|l as a 
function of ^ — ppL cos ip. Strong spin-rotation effect ex- 
pected at L > ioi produces an additional phase /? ~ 1 
which leads to splitting of the resonances as a function 
of pfL. As a result, at L > Lq oscillations discussed 
in Ref. 8 will have a twice shorter period and reduced 
amplitude. 
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V. SPECTRUM EQUATION AND CURRENT 
FOR THE JUNCTION OF AN ARBITRARY 
LENGTH 

In this section we find equations determining the An- 
dreev levels for an arbitrary length of the contact, with 
the main purpose to demonstrate, that the semiclassical 
average of the Josephson current is independent upon the 
SO coupling for any i/^o ratio. Here we use alternative 
method of calculation: instead of expression of the An- 
dreev levels via the transmission eigenvalues, we employ 
direct matching of wavefunctions obeying BdG equations 
in the 2DEG and both superconductive regions. To sim- 
plify calculations in this Section, we will consider here a 
model of equal effective masses, m„ = m. 

Eigenfunctions of the BdG equation for the S-2DEG- 
S junction can be represented as 8-dimensional vectors, 
since they contain three 2-dimensional blocks: i) elec- 
tron and hole components, ii) two spin projections, and 
iii) two direction of momentum in the x direction, px = 
±Pi?| cosiysj. Matching conditions for wavefunctions on 
both S/N interfaces consist of 16 scalar equations which 
relate 8 wavefunction components in the 2DEG region 
with 4 components in each of superconductive terminals 
(in the case of the subgap Andreev levels which decay ex- 
ponentially into the bulk of the superconductors). The 
next step is to reduce this system of 16 equations down to 
8 equations which couple together 4-1-4=8 amplitudes of 
wavefunctions in superconductors. The solvability condi- 
tion for this system of 8 linear equations is equivalent to 
the condition of vanishing of the corresponding determi- 
nant, g{e, x) = 0, which is equivalent to the one defined 
in Eq. ©. After some tedious calculations, the equation 
(/(e,x) = can be transformed into the following form: 

ff(e,x) = g+(e,x)5-(e,x) = 0, where 

5±(e,x) = cos2,f- Qcos/3± - sin/3, (26) 

where the parameter (3 is defined in Eq. H2U|) . 



= cos 5' 



Ak'^K'^/S? (cos * -f- cos x) 

(if2 _ fc2)2 _ g2) ' 

2kKe 



= 2 arctan ■ 



(^2 + /fc2)VA2 - £2 



and 



£, (27) 



with £ = 2meL/k being the energy dependent part of the 
phase ^(e). Eqs. H26() demonstrate, that in the presence of 
the Rashba interaction, the Andreev levels are generically 
spin-splitted for the contact of an arbitrary length. 

In the limit of vanishing spin-orbit interaction a = 0, 
as well as for electron trajectories with py = ip = 0, the 
spectrum equation reduces to the standard equa- 
tion cos 2^ = Q with a two-fold degenerate (due to the 
spin) solutions. In the special case of ideally transparent 
boundaries pf = Pp general spectrum equation H26() also 
reduces to the standard equation cos 2^ — Q, which now 
simplifies to 



cos + 2 arccos ^) = cos 



X- 



(28) 



For a relatively short contact with < 5 ^ 1 we 
expand the Andreev spectrum equation (|26|l in powers of 
small parameter k = "^f^ , and find the first correction 
to the result {T)) obtained in Sec. II for the k ^ limit : 



AWl 



r± sin2 1 ( 1 



3/2 



. X 

sm- 



coth X ] ,(29) 



where T± are defined in Eq. H19|l . 

In the general case of an arbitrary length of the con- 
tact, the spectral equation is too complicated to be 
solved explicitly for the energies of the Andreev levels. 
Moreover, one should remember that for a junction of an 
arbitrary L/^ ratio the continuous part of the spectrum 
(scattering states) contributes to the Josephson current 
as well as the localized levels we have considered. How- 
ever, total Josephson current (carried by both localized 
Andreev levels and continuous part of the spectrum) can 
be found following the method presented in Ref. |23j in 
terms of the spectral function g{€,x) itself. 

We will use Eqs. (1.9), (A.48) and (A.49) from Ref.lH 
modified in our case due to the presence of the spin- 
splitting, and the continuous scattering channels charac- 
terized by the transverse momentum py. Therefore the 
total current contains an integral over all py : 



W =^^Y^E / S E -9xln<?.(*-',x), (30) 



'total 



dpy 



here the summation goes over positive Matsubara fre- 
quencies: uj = 2TTT{n + 1/2), n = 0, 1, ... 

Within the semiclassical limit {Lpp ^1), calculation 
of the integral over py in Eq. H30|l can be simplified by the 
same method which was used previously in the last part 
of Sec. III. Namely, we first average over the period of the 
fast oscillations of cos^ = cos(fcL) at a fixed angle (p, 
and then do the integration over (p. Integration d^... 
in Eq. Ij^Ufl leads to the result which does not contain the 
spin-orbital phase /?: 



-ftotal(x) 



W71>0 



where Q = Q{e — iujn,x) defined in Eq. H27|l . 
and we took into account that dpy = pp coaipdip. 
Eq. (|31() demonstrates that the semiclassical average of 
the Josephson current through the S/Rashba 2DEG/S 
contact does not depend on the Rashba coupling constant, 
and this result is valid for an arbitrary Fermi velocity mis- 
match and arbitrary length of the contact. We note how- 
ever that this result is valid as long as electron-electron 
interaction in the 2DEG region is neglected. 



VI. DISCUSSION 

The above results were obtained for a model of an in- 
finitely long junction in direction perpendicular to the 



7 



current, when due to the translational invariance the 
motion along y axis was completely determined by the 
wavevector py of the corresponding plane wave. An ob- 
vious generalization of such a model system would be 
a junction with periodic boundary conditions in y di- 
rection. In this case, all our results would stay intact, 
up to replacement of continuous Py by the discrete set of 
wavevectors p„ = 2TTn/Ly. Although being somewhat ex- 
otic for SNS junctions, such a geometry does not seem to 
be impossible, taking into account recent advances in fab- 
rication of complicated InAs structures, cf. e.gm^. Usu- 
ally, however, S-2DEG-S structure is of a finite length 
in y direction Ly with closed boundary conditions, so 
that the channel eigenstates are characterized by stand- 
ing waves - mixtures of plane waves e*^"^* and e"*^*^". 
In presence of the Rashba term, the direction of the elec- 
tron momentum is coupled to the direction of its spin, 
thus determination of correct standing-wave eigenstates 
is not trivial. The major effect of a finite Ly ^ Lq will 
be the existence of a discrete set of transmission chan- 
nels, Nch = 'iLy/Xp, where Af is the Fermi wavelength 
of the 2DEG. However, some qualitative effect of closed 
boundary conditions will occur: now one would deal with 
a real basis of scattering states, thus Kramers theorem 
for the transmission eigenvaluesii would be applicable. 
It means that for closed boundary conditions and in the 
short-junction limit L/^o ~* 0; no spin-splitting of An- 
dreev levels may occur. In other terms, in a closed (in y 
direction) system Rashba coupling modifies transmission 
eigenvalues, but does not split it. How can we recon- 
cile it with a natural idea that for very long Ly the type 
of boundary conditions should not be important? The 
point is that the total Andreev spectrum of the system is 
doubly degenerate for the periodic boundary conditions 
as well as for the closed ones. In the first case, the de- 
generacy is due to the symmetry of T±{py) with respect 
to Py — > —py reflection, whereas in the second case it is 
due to Kramers theorem. In order to get global Andreev 
spectrum without degeneracy, the time-reversal symme- 
try should be broken. In particular, it happens if nonzero 
L/^o ratio is taken into account, as reported in Ref. 13 
Another possibility might be related with an open sam- 
ple geometry like the one used in Ref. Hj, where additional 
current can be passed in the direction transverse to the 
supercurrent. 

Individual Andreev levels could possible be observed 
experimentally by microwave spectroscopy or by mea- 
surement of tunnelling conductance into the 2DEG re- 
gion from an additional point-like junction. One version 
of the former type of experiment was proposed theoret- 
ically, for a single-channel junction, in Ref. figl In this 
case, resonant frequency is very high, of the order of A/?i, 
since the only possible transitions are between the posi- 
tive and the negative Andreev levels. This frequency is 
about 0.4 Thz for Nb terminals (considerably lower fre- 
quencies can be found in the case of a very small reflection 
probability, 1 — T 1 and a phase difference % « tt) . In 
many-channel junctions energy spacing between neigh- 



boring Andreev levels is reduced as 5e ~ A /Nch, but 
usually (without spin-orbit coupling) it is not possible 
to observe microwave-induced transitions between levels 
which belong to different conduction channels. The rea- 
son is the momentum conservation: different transmis- 
sion channels are characterized by different wavevectors 
Py/h which are spaced by n/Ly, whereas photon wave- 
length Xph = hc/Se is much longer than Ly, their ratio is 
of the order of {EjP^'^ / A){c/vn) ^ 10^- It seems pos- 
sible that this selection rule will not be effective in the 
considered situation with Rashba coupling, which mod- 
ifies conduction channels considerably at Ly > Lq- The 
point is that now conduction channels will be defined in 
the entangled space of orbital and spin variables, thus 
there seems to be no reason for vanishing of the inter- 
channel photon matrix element. However, this question 
will certainly need further investigation. 



VII. CONCLUSIONS 

We investigated the dependence of Josephson currents 
in clean S-Rashba 2DEG-S proximity junctions upon the 
Rashba spin-orbit interaction. We have generalized the 
Beenakker's formula for the Andreev levels for the case 
of spin-orbital scattering and found that for the case of 
an infinitely wide junction (in direction transverse to the 
current), the Andreev levels are s pin-sp litted. This re- 
sult is in agreement with papers p^all^ . where the ef- 
fect of Rashba spin-orbital interaction upon supercur- 
rent was studied either— in the case of absence of normal 
backscattering at the interfaces {pf = Pp), or in the one- 
dimensional caseiSi. We have shown that the semiclassical 
average of the Josephson current is however insensitive 
to the Rashba coupling as long as electron-electron inter- 
action in 2DEG is neglected. 

Our results show therefore that the account of the 
Rashba spin-orbital interaction for the usual model of 
SNS junction without electron-electron interaction in the 
normal region is not sufficient to explain experimentally 
observed strong suppression of the IcRn parameter with 
respect to its theoretical value. We believe that to explain 
this suppression the electron-electron interaction should 
be taken into account together with the spin-orbital ef- 
fects. We note that electron-electron interactions in both 
density-density and spin-spin channels are not weak in 
the 2DEG structures. 

A related open problem is to find average spin polariza- 
tion {Sy) in the 2DEG region, which is expected to exist 
at a nonzero supercurrent in the S-2DEG-S structure due 
to symmetry considerations, cf. e.g. Ref. 0. Note that 
a supercurrent-induced average spin polarization will in- 
duce, in the presence of electron-electron interaction, an 
effective Zeeman field which may strongly modify the An- 
dreev levels as well as the Josephson current. 
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